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For the description of an H2 molecule an effective one-electron model potential is proposed which 
is fully determined by the exact ionization potential of the £[2 molecule. In order to test the model 
potential and examine its properties it is employed to determine excitation energies, transition 
moments, and oscillator strengths in a range of the internuclear distances, 0.8 < R < 2.5 a. u. 
In addition, it is used as a description of an H2 target in calculations of the cross sections for 
photoionization and for partial excitation in collisions with singly-charged ions. The comparison of 
the results obtained with the model potential with literature data for II2 molecules yields a good 
agreement and encourages therefore an extended usage of the potential in various other applications 
or in order to consider the importance of two-electron and anisotropy effects. 

PACS numbers: 31.10.-|-z,31.15.-p,31.15.B- 



I. INTRODUCTION 

From the very beginning of quantum mechanics the 
hydrogen atom has been considered as one of the stan- 
dard model systems. The reason lies in the simplicity 
of the theoretical description of this most basic atomic 
system. On the other hand, the description of the hy- 
drogen molecule is obviously a lot more involved due to 
the much larger number of degrees of freedom. Com- 
pared to the atomic case the complexity of the molecule 
arises, e.g., from the electron-electron interaction due to 
the second electron and the anisotropy of the charge dis- 
tribution which may lead to an orientational dependence 
of a physical quantity. Additionally, there is vibrational 
and rotational motion of the nuclei and even in a Born- 
Oppenheimer approximation one has to deal with poten- 
tial curves for all electronic states and their rovibronic 
excitations. 

Consequently, it would be desirable to have a descrip- 
tion, although simplified, of the hydrogen molecule at 
hand which is of similar complexity as the one of the 
hydrogen atom. This would allow for an easy adoption 
of already existing numerical methods which were im- 
plemented for spherical one-electron problems to the de- 
scription of molecular hydrogen. But also in complex 
systems including H2 molecules like, e.g., H2 clusters or 
H2 adsorbed on surfaces a simple description of the elec- 
tronic structure is of interest. 

A second motivation becomes even more important in 
the era of fast improving computational resources which 
may make the full description of H2 molecules feasible 
even in time-dependent processes. That is, the com- 
parison of results achieved with a full calculation with 
the outcome of a simplified description of H2 which has 
atomic rather than molecular properties and accounts for 
the second electron only by screening. An analysis of the 
differences can yield the importance of the influence of 
two-electron as well as of molecular effects, like the devi- 
ation from a spherical symmetric charge distribution. 



In the context of the latter motivation a simple one- 
electron, single-centered model potential was proposed in 
a recent work [l| which deals with H2 molecules interact- 
ing with short intense laser pulses. Since the strong-field 
ionization is known to be very sensitive to the electronic 
binding energy and the exact form of the long-ranged 
Coulomb potential the proposed model potential is de- 
signed to agree in these properties with the H2 molecule. 
Thereby, the model potential can be adjusted to an arbi- 
trary internuclear distances by taking the corresponding 
value of the ionization potential. 

Regarding the first motivation, satisfying results have 
been achieved with the proposed model potential in the 
description of an H2 target in collisions with singly- 
charged ions [1]. The calculated total and differen- 
tial ionization and excitation cross sections agree well 
with literature data down to projectile velocities for 
which electron-electron effects may become important. 
Thereby, also the dependence on the internuclear dis- 
tance is examined and the nuclear motion is taken into 
account. 

The aim of the present work is to further examine the 
proposed simple single-centered, effective one-electron 
model potential and to find out why it describes the 
properties of the hydrogen molecule in the applications 
[l|, 0] to different physical systems so well. But also the 
limits of the model in the description of H2 molecules 
should be analyzed. Therefore, quantities like excitation 
energies, electronic transition moments, and oscillator 
strengths are calculated as a function of the internuclear 
distance and are compared to literature data for a full 
H2 molecule. Also, the model is used to determine pho- 
toionization cross sections and excitation cross sections 
in collisions with projectiles in order to test its applica- 
bility to different physical systems. It may be noted, that 
in the limit i? — > the model is also suitable for the de- 
scription of atomic helium, as is shortly commented on 
in the end. 

For a one-electron description of the H2 molecule also 
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other model potentials exist. To name only three, Teller 
and Sahlin [3|] discussed a two-center approach while a 
model potential for He atoms by Hartree [3] was also 
adjusted to H2 by fitting it to the correct ionization po- 
tential. It can be obtained by integrating an effective 
hydrogen atom-like charge distribution with Gauss's the- 
orem. Another widely used model is the scaled hydrogen 
atom model which treats H2 as a hydrogen atom with 
a scaled nuclear charge in order to achieve the correct 
ionization potential. The latter model is as simple as the 
one proposed in [fj] but has also the advantage that its 
wave functions are known analytically. A disadvantage 
of the scaled potential is, however, that its long-range 
behavior is not correct. It is therefore used in this work 
for a comparison of the present results with another II2 
model potential. 

The paper is organized as follows: Sec. [ll] presents the 
model potential and discusses its properties. Sec. IIIII 
considers various applications of the H2 model, namely, 
the calculation of excitation energies, transition moments 
and oscillator strengths as well as the determination of 
cross sections for photoionization and excitation in col- 
lision processes. Furthermore, the outcome of these ap- 
plications is discussed and compared to results of exper- 
iments and theoretical treatments of the full molecular 
system. Sec. IIVI concludes on the findings. Atomic units 
are used unless otherwise specified. 



expression 
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For instance, at i? = 1.4 a. u. the numerically determined 
ionization potential and /mod(Q^) given by Eq. ^ differ 
only by 0.01%. The dependence on a simplifies even 
further in the limit | a | — > where the ionization potenial 
becomes /mod (o^) — > -^h + ct and depends only linearly on 
a as can be seen in table ID 

In order to describe an H2 molecule with a fixed in- 
ternuclear distance R a certain a has to be determined 
which fulfills the requirement that /mod (a) is equal to 
the ionization potential In^ (R) of the H2 molecule at the 
considered fixed distance R. In Table U values of a which 
yield the ionization potentials of H2 for internuclear dis- 
tances R in a, range from 0.8 a. u. to 2.5 a. u. are given. 
For a fixed R the ionization potential lu^ (R) is obtained 
by subtracting the ground-state potential-energy curve of 
H2 which was very accurately calculated by Wolniewicz 
5] from the ground-state energies of Hj. Also given is 
the a value for the limit R ^ which yields the correct 
ionization potential of the helium atom 

Since the model potential can be adopted to different 
internuclear distances R with the help of a it is possible 



II. MODEL POTENTIAL 

In order to obtain a simple model for a complex sys- 
tem the right balance has to be found, i.e., a model 
which reflects the characteristics of the full description. 
It is known that, e.g., ionization processes of H2 are very 
sensitive to the ionization potential 1^2 ^-nd the prop- 
erties of the bound states depend on the exact form 
of the Coulomb potential. Hence, it is important that 
the model potential agrees in these properties with the 
molecule. An appropriate trade-off for the description 
of H2 molecules may be achieved by using the following 
simple model potential [ij] for an effective electron with 
the radial coordinate r 




where a is a dimensionless term. The model potential 
satisfies the conditions Vmodir) — > — 1/r for r 00 and 
describes therefore the long-range behavior of an effective 
H2 potential correctly as being hydrogen-atom like. Fur- 
thermore, it reduces to the potential of a hydrogen atom 
H for a with an ionization potential /h = 0.5 a. u. 

The exact dependence of the ionization potential 
/mod (a) on a for a system described by Vmod can be 
determined numerically (cf. [l|). However, an advantage 
of the model proposed in Eq. ^ is the possibility to 
approximate /mod (a) quite accurately with the analytic 
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a{R) 







0.87910 


0.903570 


0.8 


0.348416 


0.715577 


0.9 


0.302668 


0.693373 


1.0 


0.262548 


0.672753 


1.1 


0.227258 


0.653645 


1.2 


0.196111 


0.635961 


1.3 


0.168525 


0.619606 


1.4 


0.144021 


0.604492 


1.4487 


0.133081 


0.597555 


1.5 


0.122196 


0.590531 


1.6 


0.102722 


0.577647 


1.7 


0.0853182 


0.565762 


1.8 


0.0697585 


0.554815 


1.9 


0.055851 


0.544745 


2.0 


0.0434376 


0.535499 


2.1 


0.0323864 


0.527029 


2.2 


0.0225906 


0.519292 


2.3 


0.0139698 


0.512251 


2.4 


0.00646727 


0.505869 


2.5 


0.00012071 


0.500115 



Table I; Values of a used in this work for different internu- 
clear distances R in a. u. The ionization potential In^iR) 
of H2 for these 7? is also given in Hartree. It is obtained us- 
ing the H2 ground-state potential-energy curve calculated by 
Wolniewicz [5| ■ The ionization potential of a He atom [61 and 
the corresponding a value are also given as the limit _R ^ 0. 
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to study vibrational effects as was proposed in [7|, |3| ■ Ion- 
ization cross sections which account for the vibrational 
motion of the H2 nuclei in collisions of H2 targets mod- 
eled by Vmod with antiprotons were obtained in 2"! . They 
were achieved by employing closure, exploiting the lin- 
ear behavior of the ionzation cross section with R, and 
performing the calculations at R = (R) = 1.448a. u. 
(a = 0.13308). 

However, a molecule treated in the fixed- nuclei approx- 
imation differs from an atom owing to the anisotropy 
of the electronic charge distribution which cannot be 
described correctly within an isotropic, single-centered 
atomic model potential. The effect of anisotropy due to 
both the two nuclei and due to the second electron in H2 
is to some extent included as a screening of the Coulomb 
potential. Two-electron effects, like double excitation or 
double ionization, are naturally not described properly 
by the model. 

In order to compare the properties of the proposed 
model potential Vmod in Eq- © with another quite pop- 
ular (see, e.g., 0) simple artificial atomic model a scaled 
hydrogen atom Hgcai may be introduced. Its potential 

VsUr) - (3) 
r 

differs from a normal H atom due to the scaled nuclear 
charge Zscai ■ The correct ionization potential of H2 at a 
given R can be obtained for Hscai, if the nuclear charge 
is scaled as 

^scal(i?) - (/H.(i?) //h )'^' . (4) 

It may be alluded that due to the scaling of the nuclear 
charge in Eq. Q all energies ej of the bound states of 
Hscai are affected in the same way, i.e., they are shifted 
in comparison to the H atom as 

e, [Hscal] = (^scal)' P] . (5) 

Although the ionization potential of the H2 molecule is 
well described by the scaled hydrogen atom it can be ex- 
pected that this is not the case for the energies of the 
bound states, since the potential in Eq. ^ does not 
have the correct r dependence. Furthermore, one expects 
problems in the description of molecular properties that 
are very sensitive to the asymptotic long range behavior 
like tunneling ionization in intense electric or electromag- 
netic fields. 

The physical quantities studied in this work like os- 
cillator strengths, transition probabilities, or cross sec- 
tions obtained with an effective one-electron model are 
multiplied with a factor two in order to account for the 
two equivalent electrons of II2 . It should be noted, that 
also alternative ways to interprete the results of single- 
electron models for two-electron systems are possible 
(e.g. 0). 



III. APPLICATIONS OF THE MODEL 
POTENTIAL 

Since the model potential of Eq. ^ is isotropic 
it is naturally qualified for describing orientationally- 
averaged II2 molecules. This is often the case in experi- 
mental studies in which isotropic, non-aligned molecules 
are investigated. Optical excitations into p states of the 
model II2 are consequently compared with both possible 
dipole- allowed transitions into and states of the 
H2 molecule. An orientational averaging yields in this 
case the factors 1/3 and 2/3 for a weighting of the re- 
sults for the symmetries and ^n„, respectively. On 
the other hand, the isotropy is of course a limitation of 
the model. For example, in the case of multiphoton ex- 
citations interference terms prevent a determination of 
simple weighting factors [11, [Tl| . 

In what follows, it should be investigated how satis- 
fyingly the proposed model potential works in practice 
with respect to various applications. First, excitation en- 
ergies, transition moments, and oscillator strengths are 
considered. Afterwards, V,nod is used for the description 
of ionization and excitation of an II2 molecule in interac- 
tions with photons and in collisions with particles. The 
findings are compared to corresponding experimental and 
theoretical results for an H2 molecule and partly also to 
results obtained with Hgcai. 



A. Excitation energies 

In Fig. [T] excitation energies (EE) for the energetically- 
lowest dipole-allowed final states of the H2 molecule with 
the symmetries n ^I]„ and n ^Hu with n = 1, . . . , 4 are 
given in the range of internuclear distances 1 a. u. < i? < 
2.5 a. u. They were obtained from the very accurate cal- 
culations by Staszewska and Wolniewicz [H, [ll] . The 
orientationally-averaged molecular EE are given as cir- 
cles. The corresponding four EE for an atomic system 
are the energy differences Ae between the ground state 
and the 2p, 3p, 4p, and 4/ state. 

It can be seen that in all of the four cases the EE 
of the model potential approximates the orientationally- 
averaged EE of the II2 molecule very well in the whole 
R range considered here. Only for the transition into 
the 2p state the EE obtained with model potential are 
slightly higher than those for H2 for large R. It is known 
that in the R range which is considered here the 4^I]„ 
and 3^n„ states possess a dominant (ls4/) contribution 
[13, [3]. Consequently, these states cannot be compared 
to a p state of the model potential but instead to the 4/ 
state. 

In contrast to the findings for Vmod the results for 
the scaled hydrogen atom Hgcai differ from the other Ae 
curves especially for small R while they come close to 
the correct values for R > 2 a. u. For large R this trend 
could have been expected since the H2 molecule becomes 
more and more like two distant H atoms and therefore 



R [a.u.] 

Figure 1: (Color online) Excitation energies Ae of the H2 
molecule for the four energetically lowest dipole-allowed final 
states {n — 1, . . . , 4) for different internuclear distances R cal- 
culated by Staszewska and Wolniewicz p^. [T3I]: green trian- 
gles, n^'Eu states; violet squares, n^n„ states; black circles, 
orientationally-averaged '^Eu and '^IIu (see text). Present ex- 
citation energies for corresponding transitions, i.e., from the 
ground state to 2p, 3p, 4p, and 4/: red solid curve, model po- 
tential; blue dashed curve, hydrogen atom with scaled nuclear 
charge Hscai- 



Figure 2: (Color online) Electronic transition moments (TM) 
of the H2 molecule as a function of R for transitions from 
the ground state l^Eg to the four orientationally-averaged, 
energetically-lowest, dipole-allowed final states consisting of 
the symmetries and ^Hu'- black circles, Wolniewicz and 
Staszewska [l3l. [T3|: green squares, Drira Present results 
for corresponding transition moments, i.e., from the ground 
state to 2p, 3p, 4p, and 4/: red solid curve, model potential 
Vmod', blue dashed curve, hydrogen atom with scaled nuclear 
charge Ifscai. (Note the different scales, especially for the 
4/ — Is transition.) 



can be modeled by the hydrogen atom-like Hgcai- How- 
ever, it is known that, e.g., the excitation probabihty can 
depend considerably on the EE [16] and therefore should 
be described accurately, especially around the equilib- 
rium distance « 1.4 a. u. 



B. Electronic transition matrix elements 

Another test of the capability of the model potential 
Vmod given in Eq. ([1]) can be performed by considering 
transition moments (TM) which are known to be much 
more sensitive to the behavior of the wave functions than 
the energies. The dipole TM into the state | nl ) for a 
fixed i?, 

M(nO = V2 (ls|x|n/), (6) 

were computed for the same transitions which have been 
already discussed in llll Al where n and I are the principal 
and angular momentum quantum numbers, respectively. 
The factor \/2 in Eq. ([6]) accounts for the two electrons in 
the H2 molecule. TM from the H2 ground state 1 ^Sg to 
the dipole-allowed final states ^T,u and were calcu- 
lated by Wolniewicz and Staszewska [H, [l4| , Spielfiedel 
[] •") . and Drira [l3|. In Fig.[2]the orientationally-averaged 
molecular TM are compared to the present results ob- 



tained with the model potential, whereas the wrong as- 
signment done in [l3] for molecular states with dominant 
(ls4£) or (ls4/) configuration is corrected as proposed 
in O, [H, [ll|. Also given are the TM for the scaled 
hydrogen atom Hgcai- 

In general, the present TM achieved with Vmod agree 
with the data for the full H2 molecule. For R > 1.5 a. u. 
there is some deviation for the transition into the 2p 
state which could have been expected, since the electron- 
electron interaction and the effects due to the two nuclei 
are most prominent for the lowest excited states. Other- 
wise, all TM to higher states match the literature data 
very well. It may be noted that even the molecular states 
4 and 3 — both with dominant (ls4/) character 
at small R — are again nicely represented by the non- 
dipole-allowed 4/ state of the model. The EE as well as 
the vanishing TM of the 4/ state match the correspond- 
ing orientationally-averaged results of the II2 molecule. 

The TM calculated for Hgcai show for all p transitions 
a different dependence on R than the TM for H2. For 
small R they are too large and for R 2.5 a. u. they 
approach the TM calculated for Knod- The deviations 
indicate that, especially at small R, the properties of 
the wave functions obtained with Kjcai differ considerably 
from those of a real II2 molecule. 
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R [a.u.] 

Figure 3: (Color online) Oscillator strengths (OS) of the H2 
molecule as a function of R for transitions from the ground 
state l^Tjg to the four orientationally-averaged, energetically- 
lowest, dipole- allowed final states consisting of the symme- 
tries and ^n^: black circles, Wolniewicz and Staszewska 
[1^ . [13I , [1^ . Results for corresponding transition moments, 
i.e., from the ground state to 2p, 3p, 4p, and 4/: red solid 
curve, present; blue dashed curve, hydrogen atom with scaled 
nuclear charge Hgcai- 



C. Oscillator strengths 

Fig. [3] shows the oscillator strengths (OS) of the H2 
molecule, the model potential Knod, and Hgcai as a func- 
tion of R for the same transitions which were already 
considered before. It may be argued that the procedure 
of orientational averaging is most appropriate for the 
OS since they obey the Thomas-Reiche-Kuhn sum rule. 
Since the OS depend on the EE and TM the question is 
whether this leads to a compensation or even to an in- 
crease of the deviations between model and real molecule. 
The OS from the ground into the excited state | nl ) are 
given by 



f{nl) = ^ {eni - eo) \ M{nl) \ ' 



(7) 



where eo and e„; are the energies of the ground and fi- 
nal excited state \ nl), respectively. The OS for the H2 
molecule are constructed in the same way using the data 
calculated by Wolniewicz and Staszewska O [3]- 
First, the OS for both symmetries and "'^n„ were 
determined separately and afterwards orientationally 
weighted with factors 1/3 and 2/3, respectively, in or- 
der to compare to the present results. 

It can be seen in Fig. [3] that the OS of Hgcai are in- 
dependent of R and are therefore the same as for the 
hydrogen atom. This is due to a cancellation of the Zgcai 
dependence in Eq. ([7]). Therein the dependence on Zgcai 
of the energies is e oc (^scai)^ (cf- Eq. ([5])) and of the 



12 
10 



^ 8 



■ ^ 



1 1 1 


1 1 1 1 1 1 

T: present, R=1.4 - 




T: Sanchez, R= 1.4 _ 




• E: Latimer [21] 




A E: Chung [20] 




■ E: Lee [19] 


- 


^^fci^ |l 


1 


II 

1 < , , , 1 



0.75 1 

Photon energy [a.u.] 



1.25 



Figure 4: (Color online) Total photoionization cross section 
of H2 as a function of the photon energy for a fixed internu- 
clear distance R = 1.4 a. u. Theory: red solid curve, present 
results with model potential; blue dashed curve, Sanchez and 
Martin 18]. Experiment: green squares, Lee et al. JJl]; violet 
triangles, Chung et al. [20[; black circles, Latimer et al. [2l| . 



TM is M cx l/(Zscai)- It is even necessary that the OS 
of Hscai are independent of R since a scaling of the OS 
with a single factor would lead to a violation of the above 
mentioned sum rule. 

The OS of the II2 molecule and for Vnwd are, however, 
not independent of the internuclear distance R. For all 
transitions the OS of II2 and the present model agree 
well for small R. For increasing R the OS obtained with 
Knod increase roughly linearly while those for II2 show a 
different behavior for R > 1.5 a. u. However, the magni- 
tudes are still comparable. At R — 2.5 a. u. the OS ob- 
tained with Knod and Vsced coincide which already could 
have been expected before from the results for the re- 
lated EE and TM. For this distance both potentials be- 
come hydrogen-atom like. Considering the region around 
R = 1.4 a. u. — which is most important for many calcu- 
lations with fixed R considering processes starting from 
the II2 ground state — one can conclude that the OS of 
the II2 molecule are satisfyingly modeled by the proposed 
potential Vmod- 



D. Photoionization cross sections 

A calculation of the photoionization spectrum for the 
hydrogen molecule is used to demonstrate the applicabil- 
ity of the present model to interaction processes in which 
an II2 molecule is ionized. In doing so, the representation 
of the continuum states is probed. Further applications 
of Vmod in order to describe ionization of H2 in time- 
dependent processes can be found elsewhere P, Q ■ The 



photoionization cross section is given by Projectile velocity [a.U.] 

4 6 8 10 



<7ph(6) = ^(e-eo) |M(6)|%(e), (8) 

where e is the positive energy of the ionized final state 
I e ) with an angular momentum I — 1 and c is the speed 
of light. The transition matrix elements Af (e) are de- 
fined in the same way as in Eq. ^ except that the | e ) 
are considered as final states. The density of continuum 
states p(e) is used for energy-normalization of the cross 
section. 

The present photoionization cross sections were calcu- 
lated for R — 1.4 a. u. in order to compare the results 
with theoretical calculations from literature in which the 
fixed-nuclei approximation was used. Besides the theo- 
retical results by Sanchez and Martin [3] also experi- 
mental photoionization cross sections are shown in Fig.[4l 
which were measured by Lee et al. [l^ , Chung et al. |2Cll | 
and Latimer et al. [21|]. 

It can be seen that the experimental photoionization 
cross sections are well described by the present model. 
At low energies, however, the results by Chung et al. 
and Lee et al. lie slightly above and below the present 
curve, respectively. The measurements by Latimer et al. 
where performed between approximately 0.9 and 1.3 a. u. 
on a dense energy grid searching for resonances above 
1.1a. u. which they did not find. Their data match very 
well with the present curve which is, of course, free of any 
resonance caused by doubly-excited states. The clearly 
visible resonances in the theoretical data calculated by 
Sanchez and Martin around R = 1.12 and 1.25 a. u. were 
explained by Martin in (2^ as being only visible within 
the fixed-nuclei approximation. In a further calculation 
which includes nuclear motion [23| the resonance effects 
are, in accordance with experimental results, practically 
invisible. This was explained by the broadening of the 
resonances, if the nuclear degrees of freedom are included. 
For energies below 1.1 a. u. where no resonances occur in 
the data of [31 their calculation agrees well with the 
present curve. 



E. Collisional excitation cross sections 

While the ionization probability depends strongly on 
the ionization potential the excitation process is more 
sensitive to bound-state properties. Therefore, a calcu- 
lation of an excitation cross section for the H2 molecule 
is used to demonstrate that the model potential is also 
capable to describe transitions to bound states properly. 

In Fig. [5] the partial cross section for the energetically- 
lowest, dipole-allowed transition for H2 collisions with 
protons and antiprotons is shown where the H2 target is 
described with the model potential. Detailed informa- 
tion concerning the employed time-dependent method is 
given elsewhere . This transition has been chosen since 
first, it is the most probable excitation in this energy 
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Figure 5: (Color online) Differential cross sections (Jox for exci- 
tation into the lowest dipole-allowed states of H2 as a function 
of the projectile velocity v in a. u. and of the equivalent proton 
energy E in keV. Theory: Present results for a fixed internu- 
clear distance R = 1.4 a. u. for excitation into 2p. Model po- 
tential V-mod- red solid curve, antiprotons; green dash-dotted 
curve, protons. Hscai: blue dashed curve, antiprotons; violet 
dash-double-dotted curve, protons. Experiment: Sum of cross 
sections for excitations into l^Eu and l^IIu: black squares, 
electrons, Liu et al. [23I ]. 



range and second, in Sec. 1111 B] and IIII"C] the largest devi- 
ation of the TM and OS between the model and the II2 
molecule have been observed for this transition. Further- 
more, partial cross sections can be used for a more sensi- 
tive testing because the errors of total cross sections may 
be reduced by some error compensation. The present 
results are compared with experimental data for H2 col- 
lisions with electrons measured at a temperature of 10 K 
by Liu et al. [23|| due to the fact that to the best of the au- 
thors' knowledge no measurements have been performed 
for proton or antiproton impact. In addition, also results 
for proton and antiproton collisions with Hgcai are given 
in Fig. [5l 

The results for protons and antiprotons obtained with 
the proposed model potential V^nod coincide for large im- 
pact energies E > 1000 keV as is expected. At these high 
energies they also fully agree with the experimental data 
for electrons with the same impact velocity v. This be- 
havior at high impact velocities is predicted by the first 
Born approximation, i.e., the same cross section can be 
expected for collisions including particles with the same 
absolute value of the charge and the same impact veloc- 
ity. At smaller energies the cross sections start to depend 
on the properties of the projectile. Thereby, the antipro- 
ton results are closer to the measured electron data than 
the results for proton impact since the former both share 
the same charge [24] . 
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Proton energy [keV] 

Figure 6: (Color online) Electron loss cross section of He col- 
lisions with protons. One-electron loss. Theory: red solid 
curve, Vmod with a = 0.8791; blue stars, Keim et al. [2511 . Ex- 
periment: green squares. Shah and Gilbody [2^ . l27j . Net 
electron loss. Theory: red dash-dotted curve, V^od with 
a = 0.8791; blue diamonds, Keim et al. [2^. Experiment: 
black circles, Rudd et al. (28j . 



The Hgcai results for protons and antiprotons also coin- 
cide for high impact energies as expected in the first Born 
limit. However, the excitation cross sections for Hgcai in 
Fig. [5] as well as in 0] clearly show a different behavior 
than the experimental data and the results obtained with 
the potential Vmod- 

In contrast to the present results, the use of Hsced as 
a target model in calculations of total ionization cross 
sections of an H2 molecule can yield to a certain extent 
reasonable results The mixed capability of Hgcai in 
describing the H2 molecule may be explained in the fol- 
lowing way. On the one hand, concerning ionization, the 
ionization potential is in both models, i.e., Knod and 14cai, 
by definition correct. On the other hand, the potentials 
differ in their r dependence and the short-range as well as 
the long-range behavior of Vscai obviously disagrees with 
that of an H2 molecule. This leads to a poor descrip- 
tion of the bound states and finally to wrong excitation 
cross sections. This is in accordance with the deviations 
of the OS in Fig. [3] which indicate too large excitation 
probability for Hgcai at i? = 1.4 a. u. 

F. Helium atom 

In the limit i? — > with a = 0.8791 the model po- 
tential Vmod can be used for the description of a helium 
atom. Obviously, various different one-electron poten- 
tials have already been proposed in order to describe He 
atoms. There are, e.g., the well-known Thomas-Fermi 
and Hartree models as well as the Hartree-Fock-Slater 
(HFS) model [29^1 which includes a local exchange cor- 



rection and was applied, e.g., in [sO]. Another approach 
is the optimized potential method (OPM) discussed in 
(31I [3^ which was used to calculated p + He loss cross 
sections in [25'|. 

In Fig. [6] electron loss cross sections (the sum of ion- 
ization and capture) are shown for collisions of protons 
with He atoms. The present results were obtained with 
the same method which was employed for the p and p col- 
lisions with H2 in section IIIIEI and which was discussed 
in detail in [3, Measurements of the one-electron 

loss were performed by Shah and Gilbody [H, [23| . Cross 
sections for the net electron loss were experimentally de- 
termined by Rudd et al. [1^. Calculations of the one- 
electron and net electron loss were done by Keim et 
al. [2^ using the OPM with a time-independent effec- 
tive potential. 

In the present results for the net electron loss the prob- 
abilities for double capture, double ionization, and trans- 
fer ionization are counted twice in order to get the correct 
number of electrons lost in the collision process. All theo- 
retical data points for the net electron loss by Keim et al. 
coincide fully with the present findings apart from those 
for the two lowest impact energies (10 and 20keV) which 
are clearly higher than the present ones. The present 
net loss cross sections reproduce the experimental data 
by Rudd et al. to a great extent. However, in the en- 
ergy range 20 < i? < 100 keV the experimental data are 
smaller than the outcome of both theoretical investiga- 
tions. 

In the case of the one-electron loss the present findings 
match the experimental results by Shah and Gilbody well 
in the whole impact energy range of the protons. Again 
all theoretical data points by Keim et al. coincide fully 
with the present ones besides those for the two lowest 
impact energies which have again larger values. Finally, 
it can be concluded that in addition to H2 the proposed 
model potential Knpd is also capable of a simple descrip- 
tion of He atoms which is consistent with the OPM with- 
out response. 



IV. CONCLUSION 

A simple model potential Knod for an effective one- 
electron, single-centered description of the H2 molecule 
has been proposed and its properties have been exam- 
ined. The potential is unambiguously determined by the 
correct ionization potential of the H2 molecule and al- 
lows for the description of H2 at an arbitrary internu- 
clear distance R. Thereby, also the nuclear motion can 
be considered to a some extent. 

The model potential was used for various applications 
in the range 0.8 < R < 2.5 a. u. The energetically-lowest, 
dipole-allowed excitation energies, transition moments as 
well as oscillator strengths of the H2 molecule are repre- 
sented well by the present model. The model was further- 
more employed for the calculation of the photoionization 
cross section of H2 and a partial excitation cross section 
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in collisions of H2 with charged particles. In both appli- 
cations experimental and also theoretical literature data 
could be well described by the present results obtained 
with the model potential. 

Concerning the scaled hydrogen atom as model for H2 
the results for ionization are satisfying while bound states 
properties and therefore also excitation cross sections are 
not reproduced adequately. 

The satisfying description of results for H2 molecules 
justifies on the one hand the choice of the ionization po- 
tential of H2 as a criterion for adjusting Knod to a certain 
internuclear distance. On the other hand, together with 
the surpassing simplicity of the model which includes an 
approximate analytic expression for the ionization po- 
tential, it suggests its applicability to a large number of 
further problems. To name only some, there are, e.g., the 
description of the electronic structure of H2 molecules in 



H2 clusters, H2 molecules interacting with external fields 
or with particles, and finally also the modeling of He 
atoms in the limit — > 0. 

It can be concluded that the H2 molecule is in many 
cases surprisingly well described by a single-electron, one- 
center model. This means, that in these cases the effects 
of charge anisotropy and two-electron effects are small. 
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